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We develop a method for extracting accurate critical exponents from perturbation expansions 
of the 0(n)-symmetric nonlinear cr-model in _D = 2 + e dimensions. This is possible by considering 
the e-expansions in this model as strong-coupling expansions of functions of the variable e = 2(4 — 
D)/{D — 2), whose first five weak-coupling expansion coefficients of powers of e are known from 
e-expansions of critical exponents in 0(n)-symmetric (/)^-theory in D = 4 — e dimensions. 


1. Critical exponents of the 0(n)-universality class can be calculated with high accuracy from standard resummation 
procedures of renormalization group functions of (()^-field theory |^,^. For the classical Heisenberg model, where n = 3, 
the critical exponent v governing the divergence of the coherence length as ^ oc \T — Tc\~'^ has been calculated from 
seven-loop perturbation expansions in three dimensions as = 0.7073 ± 0.0030 j^], whereas five-loop expansions 
in D = 4 — e dimensions ||] extrapolated to e = 1 give v = 0.7050 ± 0.0055 Q. Apart from the initial expansion 
coefficients, the resummation procedures incorporate information on the large-order growth of the coefficients obtained 
from semiclassical considerations Results very close to the above numbers were recently obtained from a novel 

strong-coupling (()^-theory in Z? = 3 j^] as well as 4 — e dimensions |^,|^ . 

It is generally accepted that, as a consequence of the universality hypothesis of critical phenomena of all systems 
with equal Goldstone bosons, the same critical exponents should be obtainable from renormalization group studies of 
0(n)-symmetric nonlinear cr-models in = 2 -|- e dimensions at e = 1, if the second-order character of the transition 
is not destroyed by fluctuations. These conditions restrict the comparison to n > 2. For n = 1 (Ising case), there are 
no Goldstone bosons, and for n = 2 (XY-model), the transition is of infinite order, for which the divergence of the 
correlation length with temperature cannot be parametrized like ^ oc \T — , as shown by Kosterlitz and Thouless. 

Unfortunately, the e-expansions of the nonlinear cr-models have, up to now, remained rather useless for any practical 
calculation, due to their non-Borel character This has led some authors to doubt the use of such expansions 
around the lower critical dimension altogether ]. Basis of these doubts is the increasing relevance of ignored higher 
powers of the derivative term in the calculations [^ . Such a situation would be quite unfortunate, since it would 
jeopardize other interesting theories which depend on similar relationships, such as Anderson’s theory of localization 

0 - 

Fortunately, the counter-arguments are not completely convincing since they involve an interchange of limits in the 
analytic continuation in e and the increase of the number of derivatives ||l^ , so that hope remains. The purpose of 
this note is to confirm this hope and to lend further support to the intimate relationship of e- and e-expansions. This 
is done by determining from a combination of the two expansions an accurate critical exponent v for the classical 
Heisenberg model for all dimensions 2 < D < A. 


2. So far, the e-expansions of i/ ^ and the anomalous dimension 77 have been calculated up to the powers e^ 0,0: 
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The singularity at n = 2 reflects the above-discussed restriction of the upcoming considerations to n > 2. When 
evaluated at e = 1, the first series yields for the three-dimensional 0(3)-model the diverging successive values = 
(1, 2, 2.5,3.25). The often-employed Fade approximations do not help, with the best of them, the [l,2]-approximation, 
giving the too large value v = 2. So far, the only result which is not too far from the true value has been obtained 
via the Pade-Borel transform |0] 
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from which one obtains the e-dependent inverse critical exponent 

poo 

(4) 

Jo 

Its value at e = 1 is « 1.252, corresponding to ly Ki 0.799, which is still considerably larger than the accurate value 
0.705. The other Pade-Borel approximants are singular and thus of no use. See Fig. for plots of the integrands. 

A direct evaluation of the series for the other critical exponent, the anomalous dimension 77 , yields the successive 
values (2, —2, 4, —5), which are completely useless. Here the nonsingular Borel-Pade approximations [2,1], [1, 2], and 
[1,1] yield 0.147, 0.150, and 0.139, rather than the correct value 0.032. 


3. The remedy for these problems comes from a combination of the theory developed in Refs, with a procedure 
developed in Ref. j|] . The theory allows us to extract the strong-coupling properties of a (^“^-theory from perturbation 
expansions. In particular, it renders the power behavior of the renormalization constants for large bare couplings 
go, and from this all critical exponents of the system. By using the known expansion coefficients of the renormaliza¬ 
tion constants in three dimensions up to six loops, we were able to derive extremely accurate values for the critical 
exponents. The method is a systematic extension to arbitrary orders |17| of the Feynman-Kleinert variational approx¬ 
imation to path integrals p^ . For an anharmonic oscillator, this so-called variational perturbation theory yields 
expansions which converge uniformly in the coupling strength and exponentially fast, like in the order N 

of the approximation, as was observed in ]20|-|^ and proved in The extension to field theory was achieved in 

Ref. [D, and showed to same type of convergence, but with the fractional power 1/3 replaced by the irrational power 
1 — to, where u) is the critical exponent governing the approach to scaling. 

This theory is combined with the procedure of Ref. [g which allows us to interpolate variationally functions for 
which we know strong- and weak-coupling expansions. The resummation to be performed will be based on rewriting 
the above e-expansions in such a way that they may be considered as strong-coupling expansion of functions, whose 
weak-coupling expansions are provided by power series expansion in powers of e = A — D, which are known from 
(/^-theory in H = 4 — e dimensions. In this way we shall be able to derive accurate critical exponents v~^ from the 
non-Borel expansion (pT). 


4. Let us briefly recall the interpolation procedure by which a divergent weak-coupling expansion in some variable 
go of the type if at (go) = X)n=o can be combined with a strong-coupling expansion of the type Em {go) = 

9^0^ Sm=o Previously treated examples were the anharmonic oscillator with parameters p = 1/3, q = 

3 for the energy eigenvalues, and the Frohlich polaron with p = 1, q = 1 for the ground-state energy and p = 4, q = 1 
for the mass. As described in detail in [^, the first step is to rewrite the weak-coupling expansion with the help of 
an auxiliary scale parameter n as 


N 

Ajv(5o) = (^) 

n—0 

where k is eventually set equal to 1. We shall see below that the quotient p/q parametrizes the leading power behavior 
in go of the strong-coupling expansion, whereas 2/q characterizes the approach to the leading power behavior. In a 
second step we replace n by the identical expression 

K _ ^2 (0) 

containing a dummy scaling parameter K. The series (||) is then reexpanded in powers of go up to the order N, 
thereby treating as a quantity of order go- The result is most conveniently expressed in terms of dimensionless 

parameters go = go/K'^ and ct = (1 — k‘^)/go, where k = k/K. Then the replacement amounts to 

K^K{l-ago)^/\ ( 7 ) 


SO that the reexpanded series reads explicitly 


N 

W0v(5o, ct)=KP^ en{(j) {go)"", 

n—0 


( 8 ) 


with the coefficients 


2 




( 9 ) 



{p-qj)/2 

n-j 




For any fixed go, we form the first and second derivatives of WN{go, K) with respect to K, calculate the iF-values 
of the extrema and the turning points. If there is a unique extremum, this supplies us with an optimal scaling 
parameter K^. If no extremum exists, we use the turning point to determine K^. If there are more than one 
extremum or turning point, we take the smallest of these as iFjv- This procedure is called optimization. The function 
^ W^igo, Kn) constitutes the A^th variational approximation ifAr(go) to the function E{go). 

It is easy to take this approximation to the strong-coupling limit go ^ oo. For this we observe that (H) has the 
scaling form 

WNigo, K) = RP-WNigo, (10) 

For dimensional reasons, the optimal Kn increases with go like Kn w so that go = and a = 1/go = 

remain finite in the strong-coupling limit, whereas goes to zero like l/[ciq[go/■ Hence 

WN{go,KN) ~ gl^'^(?NWN{c/Yill)- ( 11 ) 


Here cn plays the role of the variational parameter to be determined by applying the optimization process described 
above to the function n {c~/Y ■, 0)- The full strong-coupling expansion is obtained from the Taylor series of WN{go, k^) 
in powers of = {go/, which yields 


WN{go) = go"^ bo{go) + higo) ^ + & 2 (ffo) 

bn{go) = ^wP{go,0)g^o'^~^^^'^ > 


with 


-4/q 


where wY {go, is the nth derivatives of WN{go, with respect to Explicitly: 
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The optimal expansion of the energy is obtained by expanding 


ffo = 7 o + 7 i 



( 14 ) 


( 15 ) 


where 70 = CjY, and finding the optimal extremum (or turning point) in the resulting polynomials 0171,72,... . 
In this way we obtain a systematic strong-coupling coupling expansion in powers of {go/i^Y This is done as 

follows: We first optimize the leading str ong -coupling coefficient bo{go) in go, and identify the optimal position by 79. 
Optimizing WN{go) with the expansion (|l^, in 71, 72,... , yields for the parameters p = — 2 , q = 2 a.t the coefficients 
71,72,... and optimal &n(do)’s by the equations listed in Table |. 

It was demonstrated in 0 how one can now find a variational perturbation series for functions for which one knows 
N weak-coupling and M strong-coupling expansion coefficients. We must merely extend the set of of coefficients 
oi,..., qn by AI unknown ones on+i, ■ ■ ■, aN+M, and determine the latter via a fit of the resulting strong-coupling 
coefficients bo, ■ ■ ■, 6m-i to the known ones. 


5. This interpolation procedure will now be applied to the perturbation expansion (|^) in 2-|-e dimensions, considering 
it as the strong-coupling expansion of a series in the variable e = 2(4 — D)/{D — 2) = 4(1 — e/2)/e = e/(l — e/2): 

Y) _ 4 T) _ 4 P ^ 

= 4:S~^ - 8-- + 16-- - 32 {[52 + 108C(3)] - [16 + 36C(3)] n + n^} - -+ ... . (16) 

n — 2 n — 2 (n — zj'^ 

The variable e plays the role of the variable go in the general formulas of the last section. The weak-coupling expansion 
of iy~^{e) in powers of e can be obtained directly from the e-expansions of Ref. |^, and has for n = 3,4, 5,1, the 
numerical form 
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n = 3 : ^-1 = 2- 0.45455 e + 0.071375 + 0.15733 - 0.52631 + ... , (17) 

n = 4 : = 2 - 0.5 e + 0.0833333 + 0.147522 _ 0.499944 + 1.47036 + ... , (18) 

n = 5: = 2 - 0.538462e+0.0955849£2 + 0.135442e3- 0.469842e^ + 1.34491e^ + ... , (19) 

n=l: = 2 - 0.333333£+0.0493827£2 + 0.158478e3- 0.539937£'‘ + 1.78954e5 + ... . (20) 


Extending these series by four more terms oq £^ + 07 E^ + as £* + 09 £®, we calculate the strong-coupling coefficients (H) 
by extremizing (0) with after identifying qq with e. The parameters (p, q) are equal to (—2, 2), as follows directly 
from a comparison of the strong-coupling powers (l -I- -I- . ..) with ([!)■ The coefficients ae , 07, ag, og are 

now determined to make bo{go), 6i(po)) ^2(50)) ^3(50) agree with (p^. The technique of doing this is described in detail 
in Ref. |^. 

In order to see how the result improves with the number M of additional terms in 0), we go through this 
procedure successively for M = 1,2,3,4. The successive additional expansion coefficients for the 0{n) univer¬ 
sality classes with n = 3, 4, 5, 1 are listed in Tables respectively. The four resulting curves for v ^(e) 

are shown in Figs, For n = 3, the successive critical exponents at £ = 1 taken from Fig. ^ are 

(i/i, 1/2, r's, JZ4) = (0.87917,0.75899,0.731431,0.712152). Their M-dependence is plotted in Fig. || as a function 
of the variable x = ® which makes them lie approximately on a smooth parabola intercepting the j^-axis at 

Voa — 0.695 ± 0.010. This extrapolated value is in good agreement with the above-quoted value « 0.705 from seven- 
loop calculations in 4 — £ dimensions The results for the other 0(4) and 0(5) universality classes are displayed 
analogously. The respective j^-values 0.735 ± 0.010 and 0.766 ± 0.010 agree well with the highest available six-loops 
results of Ref. [H, which are 0.737 and 0.767. 

As discussed above, the relation between the e- and £-expansions is expected to be restricted to n > 2, for physical 
reasons. It is instructive to see that the variational interpolation method reflects this problem at two places. First, 
the expansion coefficients in Table ^ shows a large irregularity for n = 1. Second, the successive approximations for 
v~^ in Fig. [^display no tendency of convergence with increasing order M of approximation. 

Finally, we plot our hipest (M = 4) approximations for n = 3, 4, 5 together with the large-n approximations for 
n = 00, 20, 10, 6 in Fig. M to see the change of the e-behavior for increasing n, which shows that the latter for n = 6 
is still far from the exact curve. This can also be seen in Fig. 3 of Ref. |^. 

For the critical exponent 77, the series (|^) reads in the variable e: 

,(£) = + 8"(” - - 1) [(6 - 2n^) [12 - „ - n^) C( 3 )] + ... , (21) 

whereas the weak-coupling expansion in powers of e obtained from the £-expansions of Ref. has, for N = 3,4, 5,1, 
the numerical form 


N = 3: r]/e^ = 5/242 -H 0.0183987 e - 0.0166488 -k 0.032432 + ... , (22) 

iV = 4: ?7/e^ = 1/48 -f 0.0173611 e- 0.0157657 0.029057 e^-k... , (23) 

N = 5: rj/e^ = 7/338 -k 0.0161453 e - 0.0148734 -k 0.0259628 -k ... , (24) 

N=l: ?7/e^ = 1/54 -k 0.01869?- 0.0176738 £2-^0.0386577 e^-k... . (25) 


These series can again be extended by four more terms 04?^ -k 05 £® -k oe £® -k are"^, making the strong-coupling 
coefficients bQ,bi, 62, hs, in Eq. calculated for Eqs. (p^-(|^) agree with those of the expansion of rj/e^ = 77/(2 — e)^ 
obtained from Eq. (|2l|). Here, however, we encounter problems: The 77-values from the interpolation come out too 
large by about a factor 2. Also 7 does not interpolate well. A more convenient combination of critical exponents will 
have to be found to apply this method. 
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TABLE I. Equations determining the coefficients b„(go) in the strong-coupling expansion ( [l2| ) and the associated 7 i = c„‘^Si 
in (^^ from the functions b„ = 6 ^( 70 ) and their derivatives. For brevity, we have suppressed the argument 70 in these functions. 


n 

bn 

—7™-i 

2 

62 + 71 5/ J- 571 ^0 

b\lb'l 

3 

fcs J- 7261 + 71^2 + 7172^0 + §71 ^1 + §71 

(62 J- 7161 -1- 7716 ® )/ 6 o 

4 

64 J- 7361 + 72^2 + 71^3 + (772 + 7173)60 
-^717261 + 77162 + 7717260 + |7i6® + m7i4'‘' 

(63 - 1 - 726 '/ - 1 - 7162 + 717260 ^^ 

+ 57161 ^' + 77i6o"*^)/6o 
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TABLE II. Coefficients of the successive extension of the expansion coefficients in Eq. for n = 3 determined from 
M = 1, 2, 3, 4 strong-coupling coefficients (4, 8, —16, 160) of Eq. (p^. 


n 

flg 


as 

flg 

1 

-203.827 




2 

-5.67653 

17.6165 



3 

-4.25622 

9.04109 

-15.7331 


4 

-3.80331 

6.87304 

-10.0012 

12.3552 


M 

TABLE III. Coefficients of the successive extension of the expansion coefficients in Eq. (po|) for n 
= 1, 2, 3, 4 strong-coupling coefficients (4, 0, 0, 221.096) of Eq. (|l^). 

= 4 determined from 

n 

QiQ 

ar 

as 

flg 

1 

-147.508 




2 

-7.91064 

37.1745 



3 

-4.59388 

12.3044 

-27.0837 


4 

-3.72613 

7.47851 

-12.2129 

16.9547 


TABLE IV. Coefficients of the successive 
M = 1, 2, 3, 4 strong-coupling coefficients (8, 

extension of the expansion coefficients 
-8/3, 16/3, 106.131) ofEq. (Q. 

in Eq. (poP 

for n 

= 5 determined from 

n ae 

ar 

as 


flg 

1 -108.648 

2 -10.1408 

3 -4.75598 

4 -3.57909 

60.7217 

15.1045 

7.84272 

-38.9689 

-14.1142 


21.6045 

TABLE V. Coefficients of the successive extension of the expansion coefficients 
M = 1, 2, 3, 4 strong-coupling coefficients (4, —24, 48, 3825.54) of Eq. (p^. 

in Eq. (poP 

for n 

= 1 determined from 

n aQ 

ar 

as 


ag 

1 -413.921 

2 -5.25285 

3 -442759 

4 -5.7343 

12.1104 

12450066 

13.7134 

-196950675 

-25.226 


38.0976 




FIG. 1. Integrands of the Pade-Borel transform (^) for the Bade approximants [1,3], [3,1], [2,2] and for [1,1], [2,1], [1,2] at 
e = 1. Only the last is integrable, yielding i/~^ « 1.25183 « 1/.79883. 
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FIG. 2. Inverse of the critical exponent v for the classical Heisenberg model in the 0(3)-universality class. Solid curve 
represents the interpolation result of fourth order. Lower dashed curves show iuterpolations of first, second, and third order. 
Upper short-dashed curves display, with decreasing dash length, the first three and four terms of the e-expansion (|^, respectively. 
Dotted curve is Fade Jl,2]-Borel approximations The fat dot corresponds to the seven-loop result in D = 3 dimensions, 
V — 0.7073 of Ref. 10,^. The four interpolations give {vi,V 2 tV 3 ,va) = (0.87917, 0.75899, 0.731431, 0.712152). These are 
extrapolated in Fig. g to infinite order, yielding v = 0.695. 



FIG. 3. Same ]dot as in Fig. but for the 0(4)-universality class. Fat dot represents six-loop result in D = 3 dimensions 
V = 0.737 of Ref. 0 . The four interpolations give {v]_,V 2 ,V 3 ,V 4 ) = (0.88635,0.810441,0.786099,0.768565). The extrapolation 
to infinite order shown in Fig. ^yields v = 0.735. 



e = A-D 

FIG. 4. Same plot as in Fig. ^ but for the 0(5)-universality class. There exists no Pade-Borel approxima¬ 
tion. Fat dot represents ix-loop result in D = 3 dimensions v = 0.767 of Ref. The four interpolations give 

{vi,V 2 ,V 3 ,va) = (0.89278, 0.842391, 0.820491, 0.802416). The extrapolation to infinite order shown in Fig. ^yields u = 0.766. 
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FIG. 5. Same plot as in Fig. but for the 0(l)-universality class (of the Ising modelh Again there is no Pade-Borel 
approximation. Fat dot represents seven-loop result in D = 3 dimensions v = 0.6305 of Refs. M,!!. The four interpolations give 
{i'i,V2,V3,V4) = (0.862357, 0.665451, 2.08686, 0.729231). Their failure to converge is illustrated graphically in Fig. 
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FIG. 6. The four successive approximations (i^i, t's, = (0.87917, 0.75899, 0.731431, 0.712152) for n = 3 (Heisenberg 
model) plotted as a function of a: = ® which makes them lie a smooth parabola line with the intercept v^a = 0.695 ±0.010. 

Numbers on top show extrapolated value and highest approximation (in parentheses). 
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FIG. 7. The four successive approximations {vi,V 2 ,V 3 ,Vii,) = (0.88635, 0.810441, 0.786099 0.768565) 
function of a; = M~^'^ which puts them on a smooth parabola with the intercept j/qo = 0.735 ± 0.010. 
extrapolated value and highest approximation (in parentheses). 


for n — 4 plotted as a 
Numbers on top show 











0.766316'(0.802416)' 
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0.84 
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FIG. 8. The four successive approximations {ui,U 2 ,i' 3 ,U 4 ) = (0.89278, 0.842391, 0.820491, 0.802416) for n = 5 plotted as 
a function of a: = which puts them on a smooth parabola with intercept = 0.766 ± 0.010. Numbers on top show 

extrapolated value and highest approximation (in parentheses). 



FIG. 9. The four successive approximations {vi,V 2 ,V 3 ,V 4 ) — (0.862357, 0.665451, 2.08686, 0.729231) for n = 1 plotted as a 
function of a; = M~^. They show no tendency of convergence towards the known seven-loop exponent Vaa = 0.630. 



FIG. 10. Gomparison of v~^ from the highest approximations of our interpolating resummation for the 0(n) univer¬ 
sality classes with n — 3,4,5 (counting from the top), with the values obtained from the 1/n-expansion to order 1/n^ for 
n = oo, 20,10, 6 (counting from the bottom). The n = 6 curve is still far from the exact one. 
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